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ABSTRACT 

A nonlinear Boussinesq model for the shoaling of ocean surface gravity waves on 
beaches is presented and compared to second-order finite depth theory. The spectral 
Boussinesq model of Freilich and Guza (1984) for uni-directional waves propagating 
perpendicular to a beach with straight and parallel depth contours is extended to obliquely 
propagating waves. Predictions of the shoaling evolution of a single resonant triad with 
two primary incident wave components driving a secondary wave component are 
compared to finite depth theory predictions of forced secondary waves. Results for both 
sum- and difference-interactions are presented for a range of beach slopes, incident wave 
amplitudes, frequencies and propagation directions. 

The comparisons show that there is a region (roughly between 10 and 4 m depth 
for typical swell amplitudes and frequencies) where both theories predict very similar 
growth of secondary wave components. Whereas Boussinesq theory typically predicts 
slightly smaller secondary wave amplitudes than finite depth theory, the dependence of 
the secondary wave response on the directional spreading angle of the primary waves 
predicted by both theories are in good agreement. However, pronounced discrepancies 
between Boussinesq and finite depth predictions are noted for very low (infragravity) 


frequency secondary waves on relatively steep beaches. 
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I. INTRODUCTION 

Wind-generated ocean surface gravity waves are the principal driving force for 
fluid motions and sediment transport in the littoral zone. Accurate models for the 
transformation of waves on beaches are needed to predict alongshore and rip currents, 
undertow, bedforms (e.g., bars and cusps) and beach erosion and accretion. In general, as 
waves shoal, amplitudes increase, wavelengths decrease and propagation directions 
change towards normal incidence to the beach. In addition to these linear shoaling and 
refraction effects, nonlinear wave-wave interactions cause transfers of energy to wave 
components with both higher and lower frequencies. 

Theories for nonlinear wave-wave interactions in deep water (kh>>1, where k is 
the wavenumber and h is the water depth) and intermediate depths (kh=O(1)) are usually 
based on Stokes’ (1847) perturbation expansion for weak nonlinearity (ak<<1, where a is 
the wave amplitude) (Phillips, 1960; Hasselmann, 1962; Longuet-Higgins and Stewart, 
1962). A nonlinear transfer of energy to higher frequencies occurs through sum- 
interactions: two primary wind-waves having frequencies m, and m,, and wavenumbers 
k, and k,, excite a secondary wave, with the sum-frequency @,+@, and (vector) sum 
wavenumber k, +k,. Phillips (1960) showed that at this order there are no resonances; 
secondary waves are forced waves that do not satisfy the gravity wave dispersion relation. 
Since the amplitudes are bounded, these forced waves are often called “bound waves”. 

The predicted forced secondary waves are phase locked to the primary wind 


waves, and typically distort a linear sinusoidal wave profile to shapes with peaky crests 





and broad troughs. Observations of high frequency pressure and velocity fluctuations in 
intermediate depths agree well with second-order nonlinear theory predictions (Herbers et 
al., 1992; Herbers and Guza, 1991, 1992, 1994). 


In shallow water, the Stokes-type perturbation solutions are valid only if the 





Ursell number, Ur = is small (Stokes, 1847; Ursell, 1953). While this assumption 


a 
k*h? ’ 
holds well for intermediate continental shelf depths, it usually breaks down on beaches, 
where Ur is O(1). In this shallow region (typically h<5 m), where dispersion is weak, 
nonlinearities become quite significant, and asymmetrical wave shapes develop with 


steep forward faces. 


The nonlinear shoaling region where Ur is typically O(1) is better represented by 


the Boussinesq equations (Boussinesq, 1871), which assume that nonlinearity =~ and 


dispersion (kh)*, are both weak and of the same order (i.e., Ur=O(1)). In Boussinesq 
theory, continuous nonlinear energy transfers take place between near-resonant wave 
triads. Mei and LeMéhauté (1966) and Peregrine (1967) extended the Boussinesq 
equations to a bottom sloping gently in 2 dimensions, and these equations have been used 
extensively to model the nonlinear evolution of waves shoaling on beaches. 

Freilich and Guza (1984) developed a Boussinesq model for unidirectional waves 
propagating over a bottom varying slowly in one dimension. The predicted shoaling 
evolution of wave frequency spectra agrees well with field observations (Freilich and 


Guza, 1984; Elgar and Guza, 1985, 1986; Elgar et al., 1990). Liu, Yoon, and Kirby 





(1985) developed a parabolic approximation to Peregrine’s equations, allowing for small 
directional spreading angles on a two dimensional beach, but applications to a continuous 
frequency-directional spectrum of waves have not been reported. 

Presently, no uniformly valid theory exists for shoaling waves. Whereas fully 
dispersive finite depth solutions are singular at the shoreline, weakly dispersive 
Boussinesq solutions break down in intermediate depths (kh=O(1)). Thus, both finite 
depth and Boussinesq solutions are needed to model the shoaling zone. For typical swell 
amplitudes (O(1 m)) and periods (O(10 s)), finite depth theory can be applied in the 
“outer shoaling region”, extending from roughly 20-6 m depth (Herbers et al., 1992), and 
Boussinesg theory in the “inner shoaling region”, between roughly 7 m depth and the 
outer edge of the breaker zone (e.g., Freilich and Guza, 1984; Elgar and Guza, 1985). 
Although formally the validity ranges of finite depth and Boussinesq theories do not 
overlap, for practical applications there is a region (roughly 5-8 m depth) where both 
theories are expected to be reasonably accurate. The errors in finite depth and Boussinesq 
model predictions in this transition region are poorly understood. 

Difference-interactions between two wave components (@,,k,) and (@,,k,) 
excite a third wave at the difference-frequency and wavenumber, (@, -@,,k,—k,) (for 
©, >@,), that typically causes a “set down” of the mean ocean surface under groups of 
high waves (Longuet-Higgins and Stewart, 1962). Difference-interactions are 
particularly important because they drive energetic nearshore low frequency motions 


(nominally 0.005-0.05 Hz), known as infragravity waves or “surf beat” (e.g., Munk, 








1949; Tucker, 1950; Herbers et al., 1995a). Infragravity waves cause seiches in small 
harbors (Okihiro, et al., 1993), and are believed to be the primary means by which 
sediment is suspended in the inner surf zone during storms (Holman et al., 1978; Wright 
et al., 1982; Beach and Sternberg, 1988). 

Whereas incident wind waves and the associated higher frequency secondary 
waves are usually dissipated in the surf-zone, the low-frequency, long-wavelength 
secondary waves forced by difference-interactions may reflect from shore. Longuet- 
Higgins and Stewart (1962) suggested that forced secondary waves at infragravity 
frequencies are released in the surf zone, reflect from the beach, and subsequently radiate 
seaward as free waves. Observations of the directional properties of infragravity waves 
are consistent with this mechanism (Herbers et al., 1995b), but quantitative models for 
the nonlinear transfer of energy to infragravity frequencies are presently lacking. Second 
order finite depth theory (Hasselmann, 1962, based on the Stokes expansion) predicts 
accurately the amplitudes of forced infragravity waves in intermediate depths (Herbers et 
al., 1994), but is singular near the shoreline, and thus cannot predict the amplitudes of 
free infragravity waves released in the surf zone (Herbers et al., 1995b). Boussinesq 
theory may provide quantitative estimates of free infragravity energy radiated from shore. 

This paper extends the work of Freilich and Guza (1984), by examining both sum- 
and difference-interactions between obliquely propagating shoaling waves. The 2- 
dimensional Boussinesq equations for a beach with straight and parallel depth contours 
are derived in Chapter I], and solved for a spectrum of waves in Chapter III. The 
predicted evolution of amplitudes and phases for a single wave triad are compared to 
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predictions based on Hasselmann’s (1962) second order finite depth theory in Chapter IV. 


The results are summarized in Chapter V. All Figures are contained in the Appendix. 




















II. GOVERNING EQUATIONS AND APPROXIMATIONS 
In this chapter, Boussinesq equations are derived for obliquely propagating waves 
on a gently sloping, impermeable beach with straight and parallel depth contours. 


Following Peregrine (1967), non-dimensional variables will be used as follows: 


(s,y.2)=— (xt. y%,2") t= 
0 





(u*, v*, w *) p=——— p* 


u, V, W )= 
( ) gh, pgh, 


where * indicates a dimensional variable and h, is a representative water depth. The z 
axis points vertically upwards, with z=0 corresponding to the mean surface. The x axis 
points onshore, and y points alongshore (Figure 1). The variables p, u, v, and w denote 


pressure and the velocity components in the x, y, and z directions, respectively. The free 


surface is defined by z= n(x, t) where x is the horizontal space vector (x,y). The bottom 


is defined by z = —h(x). 


The equations of motion and the continuity equation in non-dimensionalized form 


(neglecting viscous effects) are: 


Cu ou 
—+(u-V)u+ w—+ Vp = 0 2-1 
i (u-V) Ve (2-1) 





one We 
OZ 


OP is 
. = ti=0 (2-2) 


Z 


apg (2-3) 
OZ 


where V denotes the two dimensional gradient operator (22) , and u indicates the 


horizontal velocity vector (u,v). The surface and bottom boundary conditions are given 


by 
eae = 
= n at z=n(x,t) (2-4a) 
p=0 at z= (x,t) (2-4b) 
w= us at z= —h(x) (2-4c) 


Irrotational flow is assumed, which yields the relation: 


ee = Vw (2-5) 
OZ 


The continuity Equation (2-3) can be vertically integrated using the kinematic 


boundary conditions (2-4a,c) and Leibniz’s Theorem: 


7 on _ : 
V- fudz+ am (2-6) 








Following Peregrine (1967), the variables 1, p, and u are expanded in terms of 


the non-linearity parameter ¢, defined to be the ratio of wave amplitude to water depth 


e=% , 


N=en, ten, +... 
P=Pp, +€p, +6°p,t... (2-7) 


u=cU, +£°U, +... 
where p, = —Z, is the hydrostatic pressure of the undisturbed fluid. In shallow water, the 


vertical velocity component w is O(c) smaller than u, where the dispersion parameter, 


o = kh, is the product of the wavenumber and the water depth: 


w=0 ew, + €’w,+...) (2-8) 


Even waves with a large oblique incidence angle 8, relative to the shore-normal 
in deep water will refract so that the angle, 0(x), in shallow water is small and thus 


v <<u. Assuming that the waves incident from deep water propagate over a seabed with 
approximately straight and parallel depth contours, the scaling of the v component 


follows from Snell’s law: 


sin(@(x)}- ey ; 


where w and k(x)=(k(x),4) are the non-dimensionalized wave frequency and 


wavenumber vector. In shallow water this relation reduces to 





(x) = k(x)h(x)sin0, (2-9) 


Hence, the ratio between v and u in shallow water is O(kh), and we can write: 


u=eu, +€°u,+... 


; (2-10) 
V=O(Ev, +€ v) +...) 
The independent variables are scaled accordingly: 
ae ee ee Ce pans (2-11) 
Ox Ox’ OY Oy’ OZ 02’ ot ot’ 


Since the variables n;, p;,u; and w, (i=0,1,...) are O(1), the order of the terms appears 
explicitly when the relations (2-7), (2-8), (2-10) and (2-11) are substituted into the 
governing equations. 


The water depth h is assumed to be a slow function of x’. 


h = h(x) (2-12) 


where X = €x’ 1s a slow space variable, so that 


&— (2-13) 





ae . Oh* , 
Note that in dimensional coordinates, the beach slope ra 1s O(c). Thus, for typical 


Shoaling wave values of ¢ = 0.05,0 = 0.2, the beach slope is assumed to be o(1 07) , 


This is a realistic slope for a natural sandy beach. 
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First order relations for the wave-induced velocity and pressure fields are obtained 


by substituting Equations (2-7), (2-8), (2-10) and (2-11) into Equations (2-2) and (2-5) 


eu Oe) (2-14a) 
Oz' 
LAME gfe (2-14b) 
Oz’ 
Pi _ Oe) (2-14c) 
Oz’ 


It follows from Equation (2-14) that u,, v;, and p, are independent of z’. 
Equation 2-14c integrates to (using the dynamic free surface boundary condition 


(Equation (2-4b))): 
P= (2-15) 
The first order horizontal momentum equations can be obtained by substituting 


Equations (2-7), (2-8), (2-10), (2-11), and (2-15) into Equation (2-1). In component 


form: 


an ra ot = O(e) (2-1 6a) 
Ov, OM 
at Br TO) (2-16b) 











Similarly, substituting Equations (2-7), (2-10), (2,11), and (2-13) in Equation (2- 


6) yields the first order depth-integrated continuity equation: 


UTI ce toy Pe (2-17) 


Substituting Equations (2-7), (2-8), (2-10), and (2-11) in Equations (2-2) and (2-5) 


we obtain the second order relations 








2 ] 
ler AEE +O 2-18 
ee (€) (2-18a) 
CV, Ow 

2 =—!4(O(s) (2-18b) 
Oz' Oy’ 

2 ] 

2=-—14+40 : 
7 = + O(e) (2-18c) 


Vertically integrating the first order continuity equation (Equations (2-3), (2-10), 


(2-11)): 
oth eon Gas 
Ox’ = 0z' 
using Equations (2-4c), and (2-13) yields 
t fi Ou, 
WwW, (z') = —-(z ga rere) (2-19) 


Substitution of Equation (2-19) in Equation (2-18a-c) and integration with respect 
to z’ (using the dynamic free surface boundary condition (Equation 2-4b)) yields: 
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92 ? 
rao iby) [ens | = 





2-20 
a 2-208) 
2 | Ou 
v, =V, (x',y',t )- c ne Bx'By' +O(€) (2-20b) 
= (x’ ‘t')+ rhe 2 8% 10%) (2-20c) 
Pr =N2%,Y, 2 | Ox'dy’ 


where U, and V, are arbitrary velocity functions arising from the integration. 


These functions can be expressed in terms of the depth averaged second order flows U, 


and V, by vertically integrating Equations (2-20a,b): 


Loom l Cu 
u, =——|_ u,dz’=U,(x',y’,t’#—h? —— +O(e 2-21a 
2 py Litadz’UQ'y't' sh” Sop + OC) (2-21a) 
] l ju 
7 
V,=—— | v.dz’=V,(x’.y’,t'H-—h? —~— +O(e 2-21b 
ag V2e2' “VaR yt gh? SE +00) (2-21b) 


Finally, substituting the perturbation expansions (2-7), (2-8), (2-10), and (2-11), 
the first order relation (2-15) and the second-order relations (2-20a,b,c) and (2-21a,b) 


back into the horizontal momentum Equation (2-1) yields, componentwise: 





Aru 
: a.) = O(e”) (2-22a) 


| 
2 | — 
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3 
on, , ov, [2% u Se ty 2 _| 2 0¥8" (2-22) 
Ox’ 3 Ox'dy' dt’ 


and the corresponding second-order depth-integrated continuity equation is obtained by 


substituting Equations (2-7), (2-10), (2-11), and (2-13) into Equation (2-6): 


oh 
—+h——+—(nu, )+u, 2) = O(s’) (2-22c) 
OX 








II. SPECTRAL SOLUTIONS 
In Chapter II, an approximate set of governing equations was derived (2-22a,b,c) 
for waves shoaling on a beach with straight and parallel depth contours. In this chapter, 
these equations will be solved for a frequency-directional spectrum of waves. Equations 


(2-22a,b,c) can be expressed in terms of the first- and (depth-averaged) second-order 


velocity potential functions 6, and 4,: 





oy 
= ar 
u, 30%: +O(e); V, Bros 
ox’ Oy' 


Dropping the primes to simplify the notation: 


an = 
Nix +O 4x ema +Oo44 +91.94,, 1H }=O?) (3-la) 
es 1 
Nhy + Oj, +e(na + Oo + 04,94, -1 Hs] is O(e*) (3-1b) 
Nit +hd.,.. +e(n, +ho.,, +hd,,, +(n,9,,), +,,h, )= O(e*) (3-1c) 


| ateeond 


where subscripts t, x, y, X indicate derivatives. 
The lowest order wave field (n,,,) is assumed to be a linear superposition of 


nearly plane, shoreward propagating waves with frequency @, and alongshore 


wavenumber ¢ af 
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a >, pa = 85a (x)expfi(y,,(x)+ 4, »t)] 
b, = 2d a= Fie, ea EPH pa )+ y-@ »t)] 


(3-2) 


where @, = pAw and ¢, =qA¢, with Aw and Aé the separation of adjacent bands in the 


Fourier representation, and the spatial variation in phase is given by 








a +eT, (x) (3-3) 


The amplitude (a, , ) and phase (y,,) functions include a slow variation with x, owing 


to shoaling, refraction and non-linear effects. This same slow space variable X was used 
previously to scale the beach slope in Equations (2-12) and (2-13). In order for n, and 


o, to be real requires: 
Fig = Apc Vog = ~V p~a tha = af oe 
Linear terms involving 1, are eliminated from the governing Equations (3-1a,b,c) 


| _. 6 é é 6 | 
by cross differentiating y (3-1la) a (3-1c) and es (3-1b) - by (3-1c)): 
= = 1,» ? 
Diet — (hd...) = —(ho,,.). +. Dix), “a D iexxtt - (hd, ) —(n1o..).. - (hb) = Ole ) 
(3-4) 


bin = (Oia )y +4 Fan ~ Chau), + @ndiy) 2 5h raon — (6615), —Cibiedy ~ Maha), | =Ol6") 
(3-5) 











Equations (3-4) and (3-5) can be expressed as: 
oe 
M40 = hd 7 (6. a hd... +7 ix), eee ie hd yy —(n,9,, ). ws habs.) f= O(e*) 


(3-6) 


(with V the (scaled) horizontal gradient operator V = 22) ), which integrates to 


ae 7 ae 1, , 2 
Pin 7 hO ix. +5, ~ NO ox. re Oi) a D text —hd,,, -(n,4.,), a <n = F(t) + Ofe") 


(3-7) 


where F(t) is an arbitrary function of time arising from the integration which we set equal 
to zero. Substitution of the lowest order wave field (Equation (3-2) and (3-3)) into 


Equation (3-7) yields, after some algebraic manipulations: 


x x ] h. ‘ ] 3 ] h 2 
d5,,—hd>,, = a bs foe ae isa, a ed yh (3-8) 


p=-=x g=-= 


Tee y 3 2 y-mm-ndna® p &XPLi(y ace a as Jerrfl,.00+t9-0,:) row 


The left-hand side of Equation (3-8) is the wave equation for o,. To prevent 


resonant growth of d, (which would upset the perturbation expansion, Equation (2-7)) 
requires that the forcing terms on the right-hand side of Equation (3-8) do not contain any 
free oscillations. Since all these terms obey the lowest-order dispersion relation 


(Equation (3-3)), it follows that the right-hand side of Equation (3-8) must vanish. 


Ly 








Collecting like frequencies and alongshore wavenumbers and solving for the real and 
imaginary parts yields a coupled set of equations for the amplitude and phase evolution of 


the lowest order waves: 


l.. te 32]. -& 2 
Xx 4h 2@ > Bama, ae sin(Y, monn +Woa7 Wa) (3-9a) 


,=-—a 
P.qx P.q Pp 
4 h 8 eee 





. 2 2 
43 ey Ay wh, fa oh fa >». >, Beau siae COSC 5 mann +W -Vae) (3-9b) 


P.g P q p 
6 2@ 8a edi aaa 


The first term on the right hand side of Equation (3-9a) is the amplitude growth of 
the mode due to shoaling. The first and second terms on the right hand side of Equation 
(3-9b) are dispersion and refraction corrections to the wave phase. The double 
summation terms give the amplitude growth and phase changes of a mode due to resonant 
non-linear interactions of all possible triads in which the mode participates. For uni- 
directional, normally incident waves, Equations (3-9a,b) reduce to Freilich and Guza’s 
(1984) “consistent shoaling model” (Equations (16a,b) in their paper). 

Finally, we re-dimensionalize Equations (3-9a,b), and re-express them in unscaled 


6% 9? 


(physical) coordinates. Dropping the used previously to denote dimensional 


coordinates, and using Equation (3-3), we obtain 
1 oh = @ 


x 3 = ‘ 
Biss eT er ore aes a : ae a, Sin(y,. Rata a eas (3-10a) 
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IV. COMPARISON TO FINITE DEPTH THEORY 

In this Chapter the amplitude and phase evolution given by Equation (3-10a,b) are 
compared to predictions of second-order nonlinear finite depth theory (Hasselmann, 
1962). The Boussinesq equations were integrated over a plane beach for a single resonant 
triad consisting of two primary waves incident from deep water that drive a secondary 
wave component. Both sum- and difference-interactions are compared for a range of 
commonly observed beach slopes, incident wave amplitudes, frequencies, and 
propagation directions. The objective of these comparisons is to examine the similarities 
and differences of Boussinesg and finite depth theory predictions in the transition region 
from deep to shallow water, and in particular to investigate the dependence of the 
secondary wave growth on the directional spreading angle of the primary waves. While 
this approach is useful to describe the initial growth of a secondary wave component, in 
very shallow water the secondary wave amplitude grows to become comparable to the 
primary wave amplitudes and thus fully coupled spectral calculations are needed to 
predict the nonlinear shoaling evolution. 

For a single triad, we can simplify the notation used in Equation (3-10a,b) by 


introducing a single index i=1,2,3, where mode i has a frequency o"’ 


(0 <a" <@ <@®), alongshore wavenumber ¢‘” , amplitude a“, and phase y"”, and 


the triad satisfies the resonance rule: 


©? =H 4 q@® | 8 = 9) 4 g@ 
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The amplitude and phase evolution equations for this triad can be expressed (neglecting 


energy exchanges with any other components) compactly as: 


a = S$ +N sinB (4-1a) 
| | | KN | 
we =k +D°4R0- sa cosBN (4-1b) 
with: 
SY = enien (shoaling term) 
4 Oh i 
seat) 
k') = Bas (shallow water wave number) 
& 
ACS 
| Oo” . 
Be vb Ve (dispersion term) 
g 


Pee gh pi? 


(refraction term) 


20" 
3 3 w'? ‘2 ; ; ; 
NO} = 2p O_ ging (nonlinear interaction term) 
g 
Bo): 2. BO). B= yw? fe yw? = yw (biphase) 


A first-order finite difference scheme using an Arakawa A grid with a step size of 
0.1 m was used to integrate Equation (4-1). Test runs with grid spacings of 0.01 m and 
1.0 m yielded identical results, confirming the accuracy of the numerical scheme. 


Additional tests of the numerical scheme include a verification of energy conservation. 
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In second-order finite depth theory, two primary waves with frequencies wo” and 

(2) 7, (4) (2) (1) (2) : 
wo” (w? > @”’ > 0) and alongshore wavenumbers £” and ¢*” excite secondary wave 
components with the sum-frequency and wavenumber (a) + @ , 2? + 2), and the 


difference-frequency and wave number (@°?- @ , ¢ - ¢°). In finite depth theory, the 


free surface elevation is given by 


n(x.y,t)=a” cos(y (x) +2 y—@ t)+a cos(y® (x) +2 y—-0 t) 
ee a (x)a (xjeos(y (x) ty (x)+ [0 +2 f-fo [o +0 f) 
+Cal (x)a (x)cos(y (x)-y (x) + [2-2 f-[o je of) 


+other second-order terms 


(4-2) 


where the primary wave amplitudes a") a) and phases yi, yw? are given by the linear 


shoaling and refraction relations: 


ys 
| C.)c0s0..’ 
(1) (1) (i) (1) kk = |, 4. 
a\"(x) = eee = as wO(x)=ys “| (Ode ; i=1,2 (4-3) 


with subscripts 0 indicating deep water values at x = 0. C, denotes the group velocity and 
9" is given by Snell’s law. The coupling coefficients C, and C_ for the sum- and 


difference-interactions are given by (Hasselmann, 1962): 


23 








2 
(0 +o) 


ce= glk - k®| tanh(k” $k h}- (0 a 0) 





2 


OM gk Ko : KY . k® 


2 2 
a ears ee a 
2g 20 Po”? 








with k” the vector wavenumber of component 1 (k(x), £ My 

Finite depth theory predictions of primary and secondary wave amplitudes in 10 
m depth were used to initialize the Boussinesq model (Equation (4-1)). The biphase B in 
Equation (4-1) was matched to the finite depth value (0 or 2, depending on the sign of C.., 
Equations (4-2), (4-4)). 

Boussinesq and finite depth theory predictions of sum- and difference-frequency 
waves forced by a pair of swell components with frequencies 0.08 and 0.12 Hz on a slope 
0.01 beach are compared in Figures 2 and 3. The deep water incident wave amplitudes in 
these simulations are 0.2 m. In Figure 2 results are shown for two normally incident 
swell components. In Figure 3 the 0.08 Hz component is normally incident while the 
Q.12 Hz component has a deep water incidence angle of 60°. In all cases the Boussinesq 
and finite depth theory predictions show very similar secondary wave growth between 10 
and 4 m depth. Overall, the amplitudes predicted by finite depth theory are slightly 
higher than those predicted by Boussinesq theory, suggesting that finite depth theory 


overpredicts the secondary wave growth when the triad is close to resonance. Boussinesq 
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theory predictions of the biphase are generally close to the finite depth theory values of 0 
(sum-interactions) and 2 (difference-interactions). The secondary wave amplitude 
growth in the Boussinesq predictions is controlled by the biphase of the triad (Equation 
(4-12)). The growth vanishes when B is equal to the finite depth values of 0 or x. Hence, 
the Boussinesq and finite depth theories do not smoothly match in 10 m depth and the 
Boussinesq biphase appears to oscillate around, and converge to the small value which 
causes the secondary wave to grow at the same rate as predicted by finite depth theory. 

The amplitude and biphase discrepancies between Boussinesq and finite depth 
theory predictions are generally larger for difference-interactions than for sum- 
interactions (Figures 2, 3 and other examples discussed below). Errors are expected in 
both finite depth and Boussinesq theory predictions for triads involving a low 
(infragravity) frequency secondary wave component because the infragravity wavelengths 
are much longer than those of the incident swell. For example, in the model run shown in 
Figure (2a,c) the wavelength of the 0.04 Hz difference-frequency component varies 
between 250 and 150 m and thus the depth variations over a wavelength (1.5 to 2.5 m) 
are appreciable. 

While the predicted sum-frequency secondary wave amplitudes are only slightly 
smaller for obliquely propagating primary waves than for normally incident primary 
waves (cf. Figures 2b and 3b), the 60° angular separation of the primary waves causes 
about an order of magnitude reduction in the difference-frequency secondary wave 
amplitude (cf. Figures 2a and 3a). The dependence of the sum- and difference-frequency 


secondary wave amplitudes on primary wave incidence angles is illustrated in Figures 4- 
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6. Secondary wave amplitudes in 4 m depth, normalized by the product of the primary 
wave amplitudes (in the same depth), are shown versus the directional spreading angle of 
the primary waves in deep water for a range of beach slopes and primary wave 
frequencies and amplitudes. In finite depth theory the normalized amplitude is equal to 
the coupling coefficient C (Equation (4-2)) and thus can be interpreted as the secondary 
wave response. 

In Figure 4 results are shown for primary wave pairs with frequencies (0.09 Hz, 
0.11 Hz) and (0.07 Hz, 0.13 Hz), and amplitudes 0.2 m ona slope 0.01 beach. In Figure 
5 the response of a (0.08 Hz, 0.12 Hz) primary wave pair with amplitudes 0.1 m is 
compared to the response of the same frequency pair, but with larger amplitudes 0.4 m, 
on a slope 0.01 beach. Results for 0.2 m amplitude primary waves with frequencies (0.08 
Hz, 0.12 Hz) on two different beach slopes 0.003 and 0.03 are shown in Figure 6. The 
Boussinesq and finite depth theory predictions of the dependence of the secondary wave 
response on the directional spreading angle of the two primary waves are generally in 
good agreement. In all cases sum-interactions are insensitive to the spreading angle with 
approximately a 10-20% decrease in the response when the spreading angle increases 
from 0 to 60°. The magnitude of the predicted sum-frequency response is insensitive to 
the primary wave frequencies and the beach slope with Boussinesq predictions 
consistently about 25% lower than finite depth predictions (cf. Figures 4b,d; 6b,d). 
However, while Boussinesq predictions of the sum-frequency response for 0.1 m 
amplitude primary waves (Figure 5b), are in very close agreement (about 10-20% lower) 
with the finite depth values, the Boussinesq predictions for 0.4 m amplitude waves are 
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about 30-50% lower than the finite depth values (Figure 5d). In this case finite depth 
theory predicts secondary wave amplitudes that are comparable to primary wave 
amplitudes and thus is not valid. Boussinesq theory predicts a reduced response for these 
large amplitude primary waves that is virtually independent of the spreading angle. 

The difference-frequency secondary wave response is much more sensitive to the 
directional spreading angle than the sum-frequency response (cf. Figures 4a,c; 5a,c; 6a,c 
to Figures 4b,d; 5b,d; 6b,d). An increase in deep water spreading angle from 0 to 60° 
causes a reduction in the difference-frequency response by a factor 3 to 10. The 
Boussinesq predictions are generally in good agreement with finite depth theory 
predictions, showing somewhat lower response levels but approximately the same 
variation with directional spreading angle. However, the agreement is much better on a 
gentle (0.003) slope (Figure 6a) than on a steep (0.03) slope (Figure 6c). In the latter case 
the depth changes from 10 to 4 m over a distance comparable to the secondary 
wavelength and the small beach slope assumption used in both Boussinesq and finite 
depth theory is grossly violated. Similar, but less snail discrepancies, between the 
Boussinesq and finite depth theory predictions for 0.02 Hz secondary waves (Figure 4a, 
with wavelengths 0(500 m)) suggests a break down of the small beach slope 
approximation. 

Overall, the comparisons of Boussinesq and finite depth theory predictions of 
secondary wave growth in the transition region from deep to shallow water show 
surprisingly good agreement. The Boussinesq predictions of both sum- and difference- 


frequency secondary wave amplitudes are slightly lower than finite depth predictions but 
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the dependence on directional spreading angle predicted by both theories are in excellent 


agreement. 
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V. CONCLUSIONS 

Although wave transformation on beaches is generally well described by one- 
dimensional models, the directionality of waves is of crucial importance to a variety of 
nearshore processes including infragravity waves, longshore currents and sediment 
transport. In this paper Freilich and Guza’s (1984) one-dimensional spectral model for 
the shoaling of uni-directional waves on a beach, based on the Boussinesq equations for 
weakly nonlinear, weakly dispersive waves in slowly varying depth (Peregrine, 1967) is 
extended to two dimensions. A beach with straight and parallel depth contours is 
assumed on which wave incidence angles are reduced by refraction. In this 
approximation, dispersion and directional spreading effects are shown to be of the same 
order, and near-resonant nonlinear interactions occur between any pair of incident wave 
components. A two-dimensional Fourier representation in the frequency and alongshore 
wavenumber domain is used to describe the lowest order wave field, yielding a coupled 
set of equations for the amplitude and phase evolution of each Fourier mode. For uni- 
directional waves these equations reduce to Freilich and Guza’s model. 

Model results for the simple case of a single triad of waves are compared to 
second-order nonlinear, fully dispersive finite depth theory (Hasselmann, 1962) for a 
typical range of beach slopes and swell amplitudes, frequencies and propagation 
directions. In these model runs two primary swell components force a third (secondary) 
wave component with the sum- or difference-frequency and alongshore wavenumber. 
The Boussinesq model (not valid in deep water) was initialized using the finite depth 


theory predictions of the amplitudes and biphase in 10 m depth. Boussinesq and finite 
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depth theory predictions of both sum- and difference-interations are compared in the 
transition region from deep to shallow water (10-4 m depth) where both theories are 
expected to be reasonably accurate. 

In general, both theories predict very similar secondary wave amplitude growth, 
with Boussinesq predictions typically slightly lower than finite depth predictions. The 
dependence of the secondary wave response on the directional spreading angle of the 
incident swell components in deep water predicted by both theories are in excellent 
agreement. Whereas the sum-frequency response is only slightly reduced for large 
spreading angles, difference-interactions are sensitive to the primary wave directions with 
typical reductions in secondary wave amplitudes of a factor of 3 to 10 when the primary 
wave spreading angle increases from 0° to 60°. 

Not surprising, discrepancies between Boussinesq and finite depth theory are most 
pronounced for low (infragravity) frequency secondary waves on a relatively steep beach. 
In this case depth variations are appreciable over distances comparable to the secondary 


wavelength and the small beach slope assumption used in both theories is violated. 
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APPENDIX 


Figure 1. Definitional sketch of variables and coordinate frame. Obliquely 
propagating waves shoal upon a beach having straight and parallel depth contours. ho is a 
representative water depth in the region, @ is the wave propagation angle, and k is the 
wavenumber vector. 

Figure 2. Comparisons of Boussinesq and finite depth theory predictions of the 
shoaling evolution of a single wave triad. A pair of normally incident primary swell 
components with frequencies 0.08 Hz and 0.12 Hz drives a difference-frequency (0.04 
Hz) secondary wave (left panels) and a sum-frequency (0.20 Hz) secondary wave (right 
panels). The upper and lower panels show the evolution of the amplitudes and biphase. 
Deep water primary wave amplitudes are 0.2 m. 

Figure 3. Same model comparisons as shown in Figure 2 but with an obliquely 
propagating 0.12 Hz primary wave component (deep water incidence angle of 60°). 

Figure 4. Secondary wave response (defined as the secondary wave amplitude 
normalized by the product of the primary wave amplitudes) in 4 m depth en the deep 
water directional spreading angle of the primary wave components. The primary wave 
frequencies are 0.09 Hz, 0.11 Hz (upper panels) and 0.07 Hz, 0.13 Hz (lower panels). 
The lower frequency primary wave component is normally incident while the deep water 
incidence angle of the higher frequency primary wave varies between 0 and 60°. 
Difference and sum interactions are shown in the left and nght panels, respectively. The 


beach slope is 0.01 and primary wave amplitudes in deep water are 0.2 m. 
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Figure 5. Secondary wave response in 4 m depth for primary waves with 
frequencies 0.08 and 0.12 Hz on a beach with slope 0.01. Upper and lower panels show 
results for deep water primary wave amplitudes of 0.1 and 0.4 m respectively (same 
format as Figure 4). 

Figure 6. Secondary wave response in 4 m depth for 0.2 m amplitude primary 
waves with frequencies 0.08 and 0.12 Hz on beaches with slope 0.003 (upper panels) and 


0.03 (lower panels) (same format as Figure 4). 


a2 














DIFFERENCE INTERACTIONS SUM INTERACTIONS 


X = finite depth 
o = Boussinesq (a) (b) 


b- A KGHE AE SEARS SS --HOHE TS ROSES KE SF 


AMPLITUDE (m) 


7p) VEO ROBE K XA KAKA AANA 
O 


OO 
OO 
oO 140 SC0C00O0—O€ 









10 8 6 4 10 8 6 
DEPTH (m) DEPTH (m) 


Figure 2 


34 














DIFFERENCE INTERACTIONS SUM INTERACTIONS 


X = finite depth 

o = Boussinesq (a) (b) 
BOF OXO XD 3O— FO OHS KOM FF be 36 -OX-O-KO 3O- BR &K OKO HO HE —K 
BOK OX O- KOO KR: : Qe: Oy B30: SP: s 


AMPLITUDE (m) 





& 220 
= 200; .00 
O 
D 1808 ABR SOR HX SISO 
r OO 00 
oO. 140 
8 8 
DEPTH (m) DEPTH (m) 
Figure 3 


35 











DIFFERENCE INTERACTIONS SUM INTERACTIONS 


x=finite depth 
o=Boussinesq 





SECONDARY WAVE RESPONSE SECONDARY WAVE RESPONSE 





i 1 
10 
” 0 20 40 60 0 20 40 60 
PRIMARY WAVE SPREADING ANGLE (deg) 


Figure 4 


36 














DIFFERENCE INTERACTIONS SUM INTERACTIONS 





X = finite depth 
10" o = Boussinesq 





SECONDARY WAVE RESPONSE SECONDARY WAVE RESPONSE 


-1 
10 
0 20 40 60 OQ 20 40 60 
PRIMARY WAVE SPREADING ANGLE (deg) 
Figure 5 


37 





DIFFERENCE INTERACTIONS SUM INTERACTIONS 


X = finite depth 
10" oO = Boussinesq 


SECONDARY WAVE RESPONSE SECONDARY WAVE RESPONSE 





0 20 40 60 0 20 40 60 
PRIMARY WAVE SPREADING ANGLE (deg) 
Figure 6 


38 














LIST OF REFERENCES 


Beach, R., and Sternberg, R., 1988, “Suspended Sediment in the Surf Zone: Response to 
Cross-shore Infragravity Motion”, Marine Geology 80, 61-79. 


Boussinesq, J., 1871, Théorie Générale des Mouvements qui Sont Propagés dans un 
Canal Rectangulaire Horizontal, C.R. Hebd. Séane, Acad. Sci., Paris 73, 256- 
260. 


Elgar, S., Herbers, T.H.C., and Guza, R.T., 1994 “Reflection of Ocean Surface Gravity 
Waves from a Natural Beach”, J. Phys. Oceanogr., 24(7), 1503-1511. 


Elgar, S. and Guza, R.T., 1985, “Observations of Bispectra of Shoaling Surface Gravity 
Waves,” J. Fluid Mech., 161, 425-448. 


Elgar, S. and Guza, R.T., 1986, “Nonlinear Model Predictions of Bispectra of Shoaling 
Surface Gravity Waves”, J. Fluid Mech., 167, 1-18. 


Freilich, M.H., and Guza, R.T., 1984, “Nonlinear Effects on Shoaling Surface Gravity 
Waves,” Phil. Trans. R. Soc. Land, A311, 1-42. 


Hasselmann, K., 1962, “On the Nonlinear Energy Transfer in a Gravity- Wave Spectrum, 
1, General Theory,” J. Fluid Mech., 12, 481-500. 


Herbers, T.H.C., and Guza, R.T., 1991, “Wind-Wave Nonlinearity Observed at the Sea 
Floor, Part 1, Forced Wave Energy,” J. Phys. Oceanogr., 21(12), 1740-1761. 


and , 1992, “Wind-Wave Nonlinearity Observed at the Sea 
Floor, Part 2, Wavenumbers and Third-Order Statistics,” J. Phys. Oceanogr., 
22(5), 489-504. 


and , 1994, “Nonlinear Wave Interactions and High-Frequency 
Seafloor Pressure”, J. Geophys. Res., 99(c5), 10035-10048. 


Herbers, T.H.C., Lowe, R.L., and Guza, R.T., 1992, “Field Observations of Orbital 
Velocities and Pressure in Weakly Nonlinear Surface Gravity Waves,” J. Fluid 
Mech., 245, 413-435. 


Herbers, T.H.C., Elgar, S., and Guza, R.T., 1994, “Infragravity-Frequency (0.005-0.05 


Hz) Motions on the Shelf, Part I, Forced Waves,” J. Phys. Oceanogr., 24(5), 917- 
927. 


39 





Herbers, T.H.C., Elgar, S., Guza, R.T., and O’Reilly, W.C., 1995a, “Infragravity- 
Frequency (0.005-0.05 Hz) Motions on the Shelf, Part II: Free Waves”, J. Phys. 
Oceanogr., in press. 


Herbers, T.H.C., Elgar, S., and Guza, R.T., 1995b, “Generation and Propagation of Surf 
Beat”, J. Geophys. Res., submitted. 


Holman, R.A., Huntley, D.A., and Bowen, A.J., 1978, “Infragravity Waves in Storm 
Conditions,” Proc. 16th Conf. on Coastal Engineering, American Society of Civil 
Engineers, 268-284. 


Liu, P.L.F., Yoon, §.B., and Kirby, J.T., 1985, “Nonlinear Refraction-Diffraction of 
Waves in Shallow Water,” J. Fluid Mech., 153, 185-201. 


Longuet-Higgins, M.S. and Stewart, R.W., 1962, “Radiation Stress and Mass Transport 
in Gravity Waves, with Application to ‘Surf Beats’,” J. Fluid Mech., 13, 481-504. 


Mei, C.C. and LeMéhauté, B., 1966, “Note on the Equations of Long Waves over an 
Uneven Bottom,” J. Geophys. Res., 91(B7), 7343-7358. 


Munk, W.H., 1949, “Surf Beats,” Trans. Amer. Geophys. Un., 30, 849-854. 


Okihiro, M., Guza, R.T., and Seymour, R.J., “Bound Infragravity Waves”, J. Geophys. 
Res., 97, 11,453-11,469, 1992. 


Peregrine, D.H., 1967, “Long Waves on a Beach,” J. Fluid Mech., 27, Part 4, 815-827. 


Phillips, O.M., 1960, “On the Dynamics of Unsteady Gravity Waves of finite Amplitude, 
Part 1, The Elementary Interactions,” J. Fluid Mech., 9, 193-217. 


Stokes, G.G., 1847, “On the Theory of Oscillatory Waves,” Trans. Camb. Phil Soc., 8, 
441-455. 


Tucker, M.J., 1950, “Surf Beats: Sea Waves of 1 to 5 Minutes Period,” Proc. Roy. Soc., 
A202, 565-573. 


Ursell, F., 1953, “The Long-Wave Paradox in the Theory of Gravity Waves,” Proc. 
Camb. Phil. Soc., 49, 685-694. 


Wright, L.D., Guza, R.T., and Short, A.D., 1982, “Dynamics of a High-Energy 
Dissipative Surf Zone,” Mar. Geol., 45, 41-62. 


40 














INITIAL DISTRIBUTION LIST 


Defense Technical Information Center....................cccccsccscccscccscensceuccesseuceuscanccuccescencens 


Cameron Station 
Alexandria, VA 22304 - 6145 


TMU OS ae ots eras a a cn ge ckteaceeaansioca wegen de eaaisgne ss edacnnapseandemeaceameuaseeaete 


Naval Postgraduate School 
Monterey CA 93943 - 5101 


Professor T.H.C. Herbers, Code OC/Heouuuu........cccccsssssscccccsssscecccssssccccccsnsccccccsscececeacess 


Department of Oceanography 
Naval Postgraduate School 
Monterey CA 93943 - 5121 


Professor Frenzén:-Code: MAI Ei .ccsiecsssscscicdcaceaicsietitsinaistsiieli Wadeahincediavantwesavtes eateens 


Department of Mathematics 
Naval Postgraduate School 
Monterey CA 93943 - 5121 


Professor E.B. Thornton, Code OC/Tm woe. cecccccccccccssececccccccussececccseaeececesucunese 


Naval Postgraduate School 
Monterey CA 93943 - 5121 


Provessor: fC; Lippmann. © Ode OC Mh x sesacesedasscieess us snichic teat t dae cunl cutlanans cos eitssaaeeaeseas 


Department of Oceanography 
Naval Postgraduate School 
Monterey CA 93943 - 5121 


PV NAS Ns LE CON adhesin aes ic ees ends erst ate nests cshon eu da as se erecta red Hela conan 


1590 Harbor Blvd. 
Belmont, CA 94002 


4] 


